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Abstract 

For every quantized Lie algebra there exists a map from the tensor square of the algebra 
to itself, which by construction satisfies the set-theoretic Yang-Baxter equation. This map 
allows one to define an integrable discrete quantum evolution system on quadrilateral lattices, 
where local degrees of freedom (dynamical variables) take values in a tensor power of the 
quantized Lie algebra. The corresponding equations of motion admit the zero curvature 
representation. The commuting Integrals of Motion are defined in the standard way via the 
Quantum Inverse Problem Method, utilizing Baxter’s famous commuting transfer matrix 
approach. All elements of the above construction have a meaningful quasi-classical limit. 
As a result one obtains an integrable discrete Hamiltonian evolution system, where the local 
equation of motion are determined by a classical Yang-Baxter map and the action functional 
is determined by the quasi-classical asymptotics of the universal R-matrix of the underlying 
quantum algebra. In this paper we present detailed considerations of the above scheme on 
the example of the algebra U q (sl(2)) leading to discrete Liouville equations, however the 
approach is rather general and can be applied to any quantized Lie algebra. 
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1 Introduction 


The “Yang-Baxter maps” [1] are invertible maps of a Cartesian product of two identical sets X, 

U: XxX^XxX (1.1) 

satisfying the “functional” or “set-theoretic” Yang-Baxter equation [2j, 

IZ 12 0 7Z 13 o 7Z-23 = TZ-23 0 K 13 ° K 12 . (1.2) 

This equation states an equality of two different composition of the three maps TZ\ 2 , TZi 3 and 
7^23, acting on different factors in a product of three sets X x X x X (for instance, IZ 12 coincides 
with the map ( 11 . 11 ) in the first and second factors and acts as an identity map in the third one; 
the maps TZ\ 3 and 7^23 are defined similarly). 

The interest to the Yang-Baxter maps is motivated mostly by their connection to discrete 
integrable evolution equations. Important examples of the Yang-Baxter maps as well as some 
classification results were obtained in [MI]. Mention also a related consistency around a cube 
condition (12] , which is in many cases can be associated with some set-theoretic Yang-Baxter 
equation. One common limitation of the existing methods is that Hamiltonian structures of 
the arising maps generally remain unclear. As a result there are no regular procedures for 
quantization of the Yang-Baxter maps or their applications to Hamiltonian evolution systems. 

In this paper we address these problems and present a new approach to the Yang-Baxter 
maps, which is based on the theory of quantum groups mm and naturally connected to 
Baxter’s commuting transfer matrices m and the Quantum Inverse Problem Method [16]. The 
key role in our approach will be played by the universal R-matrix [13]. This notion is associated 
with the so-called, quasitriangular Hopf algebras , which include a very important class of the 
(/-deformed (affine) Lie algebras and super-algebras (more precisely, their universal enveloping 
algebras). Let A be a quasitriangular Hopf algebra, then there exists an invertible element 
R e A (x) A, belonging to the tensor product of two algebras A , called the universal R-matrix, 
which by construction satisfies the quantum Yang-Baxter equation 

Rl2 Rl3 11-23 = R-23 R-13 R-12 , R E A ® A . (1-3) 

Similarly to Eg. (11.21) the indices here indicate how the universal R-matrix is embedded into the 
triple product A ® A ® A (for example, R 12 £ A ® A ® 1 , and similarly for R 13 and R 23 ). 
Interestingly, the universal R-matrix allows one to construct a map from the tensor square of 
the algebra A to itself mm, 


7 Z : A® A^> A® A, 


(1.4) 


acting as 

7 Z : (x®y) (x'®y') = K(x®y)R~ 1 , x,ycA. (1.5) 

It will automatically satisfy Eg. (11.21) in virtue of (11.31) and commutation relations of the quantum 
algebra A. Thus, for any quasitriangular algebra A the map ([1.51) is a solution the set-theoretic 
Yang-Baxter equation (11.21) . where the algebra A serves as a set X. At this point it should be 
noted that a quasi-triangular algebra A is, of course, not just an abstract structureless set (as it 
is implicitly assumed in the setting of Eq. (|1.2D ). but has algebraic relations between its elements, 
which must be taken into account in order for the Eq. (jl.2|) to hold. Obviously, this is a slight 
generalization of the meaning of Eo. (11.2D . but this is a well needed generalization, since the 
most important applications of the Yang-Baxter maps to dynamical systems naturally require 
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to equip the set X with an additional structure for purposes of quantization or considerations 
of Hamiltonian systems. Indeed, we show that the map (11.511 defines Hamiltonian evolution 
equations for a discrete integrable quantum system in 2D (with discrete space and time) with 
an algebra of observables formed by a tensor power of the quantum algebra A. The functional 
Yang-Baxter equation (11.21) then implies that these evolution equation obey the quantum version 
of the consistency-around-a-cube condition [12]. Further, the map (11.51) possesses a discrete 
analog of the zero curvature representation. This allows one to construct a full set of mutually 
commuting integrals of motion using the standard approach of commuting transfer matrices m 
and the Quantum Inverse Problem Method m- The unitary Heisenberg evolution operator for 
the system are expressed through the matrix elements of the universal R-matrix. 

Remarkably, all steps of our new approach to the Yang-Baxter maps admit a meaningful 
quasiclassical limit. As a result for any quasitriangular algebra A one obtains a classical Yang- 
Baxter map which automatically possesses properties of a Hamiltonian map, since it preserves 
the tensor product structure of Poisson algebras arising in the quasiclassical limit of the quantum 
algebra A. The action functional of the corresponding classical discrete integrable system is 
determined by the quasi-classical asymptotics of the universal R-matrix. 

In this paper we illustrate the above ideas for the case of the algebra U q (sl(2)), which is the 
simplest example of a quasitriangular Hopf algebra. The quantum case is considered in Sect. 2 
and the classical one in Sect. 3. The considerations are completely parallel. For the reader’s 
convenience we preserve the same numeration of the corresponding subsections in both cases. 
Note the discrete integrable equations arising here for the algebra U q {sl{ 2)) are related to the 
discrete Liouville equation, various variants of which were previously considered in |19H22j . 

2 Quantum Yang-Baxter map for the algebra U,{sl{ 2)) 

2.1 Algebra U q (sl( 2)) 

In this subsection we briefly summarise the basic properties of the algebra U q (sl( 2)). A more 
detailed description of this algebra well suited to our purposes can be found in [23] . The universal 
enveloping algebra A = U q (sl( 2)) is generated by elements H,E,F, satisfying the relations, 

[H,E]=2E, [H,F\ = —2 F , [E, F] = (q - g” 1 ) (q H - q~ H ) . (2.1) 

Here we use a normalisation, where the elements E and F are multiplied by the factor (q — q~ l ) 
with respect to the usual choice. Below, it will be convenient to also use the element K = q H , 
for which the commutation relations (12.21) become 

KE=q 2 EK ;, K F = q~ 2 F K , [E, F] = (q — <? _1 ) (K — AT -1 ) . (2.2) 

The quadratic Casimir operator, which commutes with all other elements of the algebra, has 
the form 

C = q- 1 K + qK~ l + EF . (2.3) 

The algebra (|2.2I) is a Hopf algebra with the co-multiplication A, the co-unit e, and the 
antipode S. The co-multiplication is a map from the algebra A to its tensor square 

A: A^A®A : (2.4) 

defined as 

A {K) = K®K, A(E) = E® K + \®E , A(F) = F (g> 1 + A' 1 ® F . (2.5) 
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The co-unit e is a map from A to complex numbers. It is defined as 


e(K) = l, e(E) = e(F) = 0, (2.6) 

The co-multiplication and co-unit define algebra homomorphisms, i.e, A (ab) = A (a) A (b) and 
e(ab) = e(a)e(b). The antipode S is defined atQ 

S(K) = K~\ S(E ) = —E K~ 1 , S(F) = -K F. (2.7) 

Correspondingly, for the element H one has 

A(H) = 77 01 + 10R, e{H) = 0, S(H) =-H, K = q H . (2.8) 

Note that the antipode is an algebra anti-homomorphism, i.e., S(ab) = S(b ) S(a). It satisfies an 
important properties 

5(1) = 1, eo 5 = e, (S®S)oA = aoAoS = A'oS (2.9) 

where A' is another co-multiplication obtained from A by interchanging factors in the tensor 
product, 

A' = a o A, a(x 0 y) = (y 0 x). (2.10) 


2.2 Universal R-matrix 

The algebra JJ q (sl( 2)) is a quasi-triangular Hopf algebra. This means that there exists an element 
R e A 0 A, called the universal R-matrix, which satisfies the properties 

A'(x) R = R A(x), V xe A, 

(A®1)R = R 13 R 23 , (2-11) 

(10A)R = R 13 R 12 , 

where R 12 = R 0 1, R 23 = 1 0 R and R 13 = (a 0 1) R 23 . Together with (12.71) and (12.91) this 
definition implies the quantum Yang-Baxter equation 

Rl2 Rl3 R23 = R23 Rl3 Rl2 5 

and two simpler relations 

(e 0 1) R = (1 0 e) R = (101), 

(5 0 1) R = (1 0 S -1 ) R = R 1 . 

Note, that the universal R-matrix is not unique as an element of A® A. However, if one assumes, 
that R 6 B _ 0 B + , where B + and B _ are Borel subalgebras of A, generated by the elements 
(.H , F) and ( H , E) respectively, then the universal R-matrix is uniquely defined by the following 
formal series [13] 

00 

R = q(l~q 2k+1 E®F^J . (2.14) 

k =0 V 

1 More generally, the defining property of the antipode reads [23], 

m o (5 (x) 1) o A(a) = m o (1 ® 5 _1 ) o A(a) = 0 

where m is multiplication map m o (a ® b) = ab for algebra A. For A = U q (sl( 2)) this implies (12.711 . 


( 2 . 12 ) 


(2.13) 
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2.3 Quantum Yang-Baxter map 

Let X = {K,E,F} denotes the set of generating elements of the algebra A and X\ 2 denote 
these sets in the corresponding components of the tensor product A ® A, 

X x = X®l, X 2 = 1®X, X = {K,E,F}. (2.15) 

Let us now explicitly calculate the map (11.511 . Evidently, it is completely determined by its 
action on the generating elements of the tensor product, 

U : {X 1 ,X 2 )^ (X[,X' 2 ), X[ = RX, R 1 , * = 1,2. (2.16) 

Using (|2.14p and (12.21) one reproduces the result of [18J 

' K[ = R\ (1 - q - 1 Kx~ x ExF 2 K 2 ), 

< E[ = E x K 2 , (2.17a) 

t F[ = F\ K 2 ~ 1 + F 2 - K x ~ 2 F 2 (1 -qK^ExF^y 1 , 

and 

' K’ 2 = (l- q - 1 K^ExF^y 1 K 2 , 

< E 2 = K\ E 2 + Ex- Ex K 2 2 (l - q Kx~ l ExF 2 K 2 )~ l , (2.17b) 

. K = Rx 1 f 2 . 

Note, that the relatively complicated expressions for F[ and E 2 above can be easily obtained by 
combining the remaining four equations in (12.171) with the formulae for the co-multiplications 
A and A 7 . 

For further references present also the inverse map. Solving ()2.171) with respect to the 
“unprimed” variables, one obtains 

r Kx — K[ (l — q^ 1 E[ F^j 1 , 

< Ex = E’yK' 2 y 1 (l^q- l E’ 1 F! 2 y\ (2.18a) 

. Fx = (F[ + (K'x)- 1 F 2 ) (1 - q” 1 E[ F') I<’ 2 - I<[ F 2 K' 2 , 

and 

'K 2 = (1 - q - 1 E'x F 2 ) K' 2 , 

^ E 2 = (E' 2 +E l xK' 2 ){l-q- 1 E' l F! 2 ){K' i y l -E[{K l l K' 2 y\ (2.18b) 

. F 2 = K[F'(l-q- 1 E[F'y 1 . 


2.4 Properties of the quantum Yang-Baxter map 

It is instructive to reformulate the properties of the universal R-matrix, stated above in (12.111) 
(12.131) . as properties of the quantum map ([2 nfl For this purpose is it useful to introduce a 


2 As noted above this map was previously obtained in [Tj3, however, its properties, given below in Eos. 112.281) 

(12.321) are new. 
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’set-theoretic multiplication” 6 which acts on two sets of generating elements 


<5: (X 1 ,X 2 )~X' } 


and write it as 

X' = 5(X 1 ,X 2 ). 

Explicitly, it is defined as 


K' = K\ K'2 , E’ = E 1 K 2 + E 2 , F' = F 1 + K x ~ x F 2 , 


(2.19) 

( 2 . 20 ) 


( 2 . 21 ) 


which are essentially the same formulae as in (12.51) . but their meaning is rather different. Indeed, 
it is worth noting that the direction of the arrow in (12.191) is reversed with respect to (12.41) and 
that the map (12.191) is only defined on two sets of generating elements. 

Below, it will also be convenient to write the Yang-Baxter map ()2.16[) in a functional form 

(X[,X! 2 ) = K(X 1 ,X 2 ). (2.22) 


Fig. Q] shows a graphical representation of the maps (|2.20[) and (|2.22l) . Note that these map are 
not symmetric upon exchanging their arguments. To avoid an ambiguity their first arguments 
on the diagram are marked with heavy dots. 




(X[,X' 2 ) = K(X 1 ,X 2 ) 


Figure 1: Graphical representation of the maps (12.201) and (12.221) . The first arguments are marked 
with with heavy dots. 

Consider the tensor product .40.4.0.4. of three algebras A and let (X\, X 2 , X 3 ) denote three 
sets of generators in the corresponding components of the product. Define a functional operator 
77-12 

KuiX^X^Xz) := (X[,X’ 2 ,X 3 ) = (V,(X U X 2 ),X 3 ) (2.23a) 

which acts as (|2.16l) on the first two sets and does not affect the third one. Similarly, define 
operators 

77-13 = T23 O fci2 o (723, 77 23 = <D 2 O 77 13 O <7 12 . (2.23b) 

At this point it is worth commenting on a relationship between the operator notations used in 
Sect l2.2l and the set-theoretic substitutions TZjj defined here. Note, that successive similarity 
transformations acting on a list of generator sets (X \, X 2 , Y 3 ...) with quantum operators Ry 
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is equivalent to a successive application of the corresponding set-theoretic map 7 Zij taken in the 
reverse order. For instance, one can easily verify 


r 12 r 13 (X 1 ,X 2 ,X 3 ) r^ r£ = nu(Ri2(x 1 ,x 2 ,x 3 )), 

(2.24) 

-^12-^13-^23 (Yi,X 2 , Y 3 )-R 23 i2 13 = lZ 23 (1Zl 3 (1Zl 2 (Xi,X 2 ,X 3 ))) . 

Ea. (|2.12l) then immediately implies that the set-theoretic maps 7 Zij satisfy the functional Yang- 


Baxter equation 

7^23(72-13(7^12 (Xi,X 2 ,X 3 ))) = 1Zi 2 (1Zi 3 (TZ 23 (Xi, X 2 , X 3 ))) . (2.25) 

and, more generally, 

7^23 0 Rl3 0 7^12 = 72-12 o 72-13 0 72-23 • (2.26) 

To reformulate the remaining properties of the universal R-matrix introduce additional no¬ 
tations, 

M^i,Y 2 ,Y 3 ) := (5(X\, X 2 ), X 3 ) , 5 23 (X 1 ,X 2 ,X 3 ) := (X u S(X 2 , X 3 )) . (2.27) 

Then Eas. (|2.11l) lead to 

5(X 2l X 1 ) = 5o'R(X 1 ,X 2 ), (2.28) 

72(<5i 2 (Y 1 ,Y 2 ,Y 3 )) = <5i 2 (72 23 (72i3(Y 1 ,Y 2 ,Y 3 ))), (2.29) 

72((5 23 (Xi,X 2 ,X 3 )) = 5 23 (R 12 (1Z 13 (X U X 2 ,X 3 ))) . (2.30) 


Y 2 


X’l 


X 2 




Figure 2: Graphical representation of the relations (12.2811 and (12.301) . 















Further, from (12.131) it follows that 


IZo (e®l) = (e® 1 ), 

IZo (l®e) = (1 ® e), 
Ko(S®S) = (s®s)on- L 


(2.31) 

(2.32) 


where the co-unit e and the antipode S are defined in (12.71) . We would like to stress that the 
above relations (|2.28I) - (I2.32I) are not specific to the algebra U q (sl( 2)) and must hold for any 
quasi-triangular Hopf algebra. In the case of U q (sl( 2)) they could be verified using the explicit 
form of the map TZ and its inverse given in (12.171) and (12.181) . The equations (12.281) and ()2.301) 
are illustrated in Fig. [2] 

2.5 R-matrix form of the U q {sl{ 2)) defining relations 

As remarked before the universal R-matrix is not uniquely defined by the relations (12.111) . One 
can easily to check, that if R 12 £ £>- ® £>+, given by ()2.14p . satisfies (12.lip , then so does the 
element 

R i2 = R 21 1 e B+ (x) B- , (2.33) 

which does not coincide with R 12 (this is why the algebra U q (sl( 2)) is called a quasi-triangular 
Hopf algebra). In particular, the new element R* satisfies the Yang-Baxter equation 

R 12 R 13 R 23 = R 23 R *3 R 12 j (2.34a) 

which is a simple corollary of (12.121) and (12.33[) . Similarly, one can derive “mixed” relations 


R 12 R 13 R 23 

= R 23 R 13 R 12 j 

(2.34b) 

R 12 R 13 R 23 

= R 23 R 13 R 12 , 

(2.34c) 

n d* D* 

Rl 2 Rl3 R23 

n* n* n 

- K 23 R 13 Rl2 , 

(2.34d) 

T3 * T> * n 

r^l2 R13 

n to * to * 

- K -23 R 13 14 12 • 

(2.34e) 


Let 77 1 denotes the 2-dimensional (spin ^) representations of the algebra (12.11) . 


Y S >-( l -l) ' = («-«“') (2 J). L 2 


Define the L-operators 


L + = (t7 1/2 ® l) (R) , L = {t7 1/2 ® l) (R*) , 


(2.35) 


(2.36) 


by evaluating the universal R-matrix in the 2-dimensional representation in the first space (to 
be called an auxiliary space). Using (12.141) and (12.331) . one gets 


L+ = 


H_ H_ 

q 2 q 2 F 


L“ = 


0 q 


_H_ 

q 2 0 

^ _H_ H 

—Eq 2 q: 


(2.37) 
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which are operator-valued matrices, whose elements belong to U q (sl{2)) and act in the “quantum 
space”. Choosing the two-dimensional representation for the latter 


R+ = (1 ®tt 1/2 )(L + ) = (vr 1/2 <8>7T 1/2 )(R), R = (1®7T 1/2 )(L ) = (tt 1/2 ® tt 1/2 )(R*) , (2.38) 
one gets the block i?-matrices 

/ q i \ ( (f 1 i \ 


R + = q~ 2 


1 


V 


(v-q- 1 ) 

1 


, R~ = (? + 2 


q 


i 


('q -q) 


\ 


q- 1 ) 


(2.39) 

where the internal blocks act in the second space of the product C 2 ® C 2 . It is convenient to 
define the parameter-dependent L-operator 


L(A) = AL + - A -1 , (2.40) 

and the R-matrix 

R(X) = \R + - A -1 R~ . (2.41) 

The latter, of course, coincides with the R-matrix of the six-vertex model [24] . It possesses an 
important property 

R(q*) = (q - q- 1 ) P, (2-42) 

where P is the permutation matrix, interchanging the factors of the product C 2 ®C 2 . It is useful 
to define 

R± = R± p, R(A) = R(A) P . (2.43) 

Let us now specialize the Yang-Baxter equations (12.121) . and (I2.34p by choosing the 2-dimensional 
representations (12.35[) in the first and second spaces and then using (12.361) and (|2.38p . The matrix 
elements of the L-operators will then act in the third space, which now becomes the quantum 
space. As a result one obtains the “R-matrix form” [25] of the commutation relations (12.11) 

R± (L+ <g> L+ ) = (L+ <g> L+ ) R± , (L" ® L“ ) = ( L" ® L- ) R± , 

(2.44) 

R+ (L+® L“) = (L-® L+)i?+, R~ (L^® L+) = (L+® L") Rf 2 . 

The symbol ® here denotes the tensor product of two-by-two matrices L-, acting in different 
factors of C 2 ® C 2 , and operator product for their matrix elements. With an account of (I2.42p 
and the definitions (12.401) . ([2.4ip . the above relations can be combined into a compact form of 
the parameter-dependent Yang-Baxter equation 

R( A/m) (L(A) ® L(/i)) = (L(/x) ® L(A)) R( A//x), (2.45) 

which is well-known in the theory of the six-vertex model (see [261127] for recent advances in this 
area). Finally, note that the co-multiplication A, defined in (12.51) . can be written as 

A(L±) = L± ® L ± (2.46) 


where the symbol ® denotes the matrix product of the two-by-two matrices (12.371) and the 
tensor product for their operator-valued elements, belonging to different copies of the algebra 
U q (sl(2)). There are two equations in ([2.461) . where the superscripts are either all pluses or all 
minuses, simultaneously. 
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2.6 Heisenberg-Weyl realisation 

The algebra U q (sl(2)) admits an important homomorphism into the Heisenberg-Weyl algebra, 
W q , generated by two elements u and v, 

W q : uv = ? 2 vu. (2.47) 

Namely, if one sets 

K = u, E = v (z — q u), F = v” 1 (1 — gz -1 u~ 1 ), (2.48) 

where the element z commutes with u and v, then all commutation relations (12.211 will be 
satisfied in virtue of (12.471) . The commuting element z parametrises the Casimir operator (12.311 . 

C = z + z- 1 . (2.49) 


Let us now consider the representation (|2.48l) as a change of variable for the set of generators 
{K. E, F} to the set {u, v, z} and rewrite the map (12.171) in the new variables. Assume the same 
meaning to the subscripts 1,2 as in (12.1711 (they distinguish components of the tensor product 
.4.0.4,), for instance, 


u) = 77(ui) = R (u (x) 1) R 1 , U 2 = 77 ( 112 ) = R (1 0 u) R 1 , 
and similarly for V ].2 and z\ )2 - From (12.171) it follows then 


77(ui) = uig, , 
77(vi) = (v- 1 + (v^ 


77(u 2 ) = g 1 u 2 


1 —— v 2 X ) u 2 1 
Z2 


77(v 2 ) = —v x + (v 2 —— v x ) u x , (2.51) 

Z2 Z2 

n{z 2 ) = z 2 . 


q 


77(^i) = z 1 , 
where we have used the notation 

g = 1 - (z x - qu x ) v^ 1 (u 2 - qz^ 1 ) . 


(2.50) 


(2.52) 


Note that the commuting elements z\ and z 2 , obviously, remain invariant. They enter the map 
(12.5111 as (spectral) parameters. Similarly for the inverse map, one obtains 

r 77“ 1 (ui) = mg- 1 , 77- 1 (u 2 ) = gu 2 , 

< 77 _1 (vi) = (— v^ 1 + (vf 1 - — vj 1 )u 2 > ) , 77- 1 (v 2 ) = v x + (v 2 - — v^u^ 1 , (2.53) 

V z 2 q / q 

77- 1 (zi) = zi, 77 _1 (z 2 ) = %2 ■ 

where 

g = 1 - ? -1 v 1 (zi - guj) v^ 1 (1 - qz^ui 1 ) . (2.54) 
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Q i 



Q 2 


Figure 3: (a) graphical illustration of the map 77, defined by (12.55|) . showing the assignment of 
the algebras Ai,A 2 , A [, A' 2 to the edges of a quadrilateral, (b) graphical illustration for the zero 
curvature relation (12.67bl) . The L-operators play a role of edge transition matrices for an auxiliary 
linear problem. 

2.7 Discrete quantum evolution system 

Here we introduce a (1 + l)-dimensional integrable quantum evolution system with discrete 
space and time. In doing this we follow the scheme suggested in [28LI29] and further developed 
in |30H32j . Remind that the quantum Yang-Baxter map (12.16(1 acts on a tensor square of the 
algebra A, 

77: (Ai®A 2 ) >-► (A' 1 ®A , 2 ) = R(A 1 ®A 2 )B .~ 1 . (2.55) 

It can be represented graphically as in FigGKa). We use this map to define a discrete quantum 
evolution system, where the algebra of observables 

O = A\ ® A 2 <8> • • • <8> A 2 N -1 <8> A 2 N 1 N>1, (2.56) 

is formed by a tensor power of the algebra A (all factors here are just independent copies of the 
algebra A; the indices indicate positions of the corresponding factors in the product). Let 

77 ij = Gij o 77 ij (2.57) 

denotes a composition of the map (|2.55|) . acting on the i-th and j-th factors of the product 
(12.56p . followed by the permutation operator , swapping these factors with each other. With 
this notation we can define the map 

U = So [tZ \2 o 77-34 o ■ ■ ■ o TZ{ 2 N- 1 ), 2 n) , (2-58) 

where the operator S cyclically shifts the factors in the product (|2.56l) . 

S : An —> An+ 1 , A 2 n+i = Ai. (2.59) 

The map (12.581) defines an evolution for one step of the discrete time. By construction it 
preserves the tensor product structure of the algebra of observables 

O = 7/(0) = Ai ® A 2 ® ® ^2^—1 ® “^ 2 /v • (2.60) 

where each algebra A' n is isomorphic to the algebra A. So the evolution map is an automorphism 
of the algebra of observables, which preserves its tensor product structure (often referred to as 
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an “ultra-local” structure). For a geometric interpretation consider the square lattice drawn 
diagonally as in Fig. 01 with 2N sites per row and periodic boundary conditions in the horizontal 
(spatial) direction. Now assign the factors of the product (|2.56l) to the set of connected edges 
forming a horizontal “saw”, as shown in Fig. [5J In the same way assign the algebras A' n in the 
product of (12.601) to edges of the saw, shifted by one time step above from the initial saw. Note 




time 

<— (t + 1 ) 


Figure 4: Assignment of dynamical variables to edges of the square lattice. Due to periodic 
boundary conditions the points Pl , P %, P-j on the left side should be identified with their images 
P[, P 2 , P 3 on the right side of the lattice. 

that these conventions are completely consistent with the graphical representation of the map 
(I2.55p . corresponding to one quadrilateral, shown in Figj3))a). 

The elements of the algebra p.56l) constitute the set of dynamical variables of the system at 
any fixed value of time. With the Heisenberg-Weyl realization (|2.48l) the generating elements of 
this algebra are expressed through the elements {uj, v,, Zij, where the index i = 1 , 2 ,... 2 N is 
considered as the spatial coordinate, numerating the factors in (|2.56l) . Combining (12.511) with 
(I2.58P and (I2.60p . it is easy to see that 


U(Z2n+1 ) " 2 n —1 j Ui,Z2n) Z2ni Tl lj 2, . . . , N , (2.61) 

where Z 2 N +1 — z i ■ For simplicity, consider the homogeneous case, when 

Z2n-l=Zl, Z 2 n — z 2 i Vn, (2.62) 


with arbitrary zi t 2 - Then all z’s stay unchanged under the time evolution, while the remaining 
variables transforms as 


U (U2n+l) 

^2n—1 Sn i 

U(y 2 n+l) 

= (Vn + ( v : 

^(u 2n ) 

= gn lu 2n , 

U(v 2n ) 

Zl 

~ - V 2n-1 + 


Z 2 


n ) ^ 2 n 


(2.63) 


Z2 


' v 2n—1) u 2n—1 J 
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where n = 1, 2,... N, 


gn = 1 -q 1 u 2n 1 _iV 2n _i ( z 1 - qu 2n -i) v 2 n (u 2n - qz 2 X ) . (2.64) 

and the cyclic boundary conditions u 2 tv + i = ui and v 2 at + i = vi are implied. With this 
realization the algebra of observables at any fixed moment of the discrete time reduces to the 
(spatially localized) tensor product of the Heisenberg-Weyl algebras (12.47|) , exactly as one would 
expect for the equal-time commutation relations of the canonical variables in a discrete quantum 
system. 


2.8 Zero curvature representation 

Consider again the Yang-Baxter map (12.551 . Let and h 2 denote L-operators (12.371 . with 
elements belonging respectively to the first and second algebra of the tensor product A\ <S>A 2 , 
appearing in (|2.55p . Next, recall the Yang-Baxter equation for the universal R-matrix, given 
in (12.121 . Cyclically shifting its indices (1,2,3) —> (3,1,2) and then choosing the 2-dimensional 
representation (12.351 in the (new) space 3, one obtains 

L+ L+ = Ri 2 L+ L+ R^ 1 . (2.65) 


This equation involves a matrix product of L’s (as two-by-two matrices) and a tensor product 
for their operator-valued matrix elements. Without the use of the indices 1, 2, the same equation 
can be written as 

L+ <g> L+ = R (1® L+) (L + <g> 1) R' 1 , (2.66) 

however we will prefer the indexed form (|2.65p for greater clarity. Further, let and are 
defined by the same equations 02.37[) , but with elements belonging, respectively, to the algebras 
A!\ and A 2 , appearing in (|2.55p . Then combining ()2.651 and (12.551 . one arrives to the equation 


L+ L+ = L+ L+ , 

Similarly, from (12. 34b j) and (12. 34c j) one obtains 

^2 = ^2 , L 1 L 2 = L 2 L x . 

For instance, the hrst equation on the last line reads 


(Ki K 2 )~l 


' K 2 k 2 f 2 

-q E x K 2 K x (1 - q - 1 K f 1 K 2 E 1 F 2 ' 


= (K[ 


k’ 2 ) 


1 — q~ l E[ F 2 ) Kg K\ K! 2 F! 2 

—q E'i 


K[ 


(2.67a) 


(2.67b) 


( 2 . 68 ) 


This equation, as well as the other two matrix equations in (12.671 are simple corollaries of the 
explicit expression of the Yang-Baxter map, given in (12.171 . Conversely, the equations (|2.671 
can be used as an alternative definition of the Yang-Baxter map independent of the notion of 
the universal R-matrix. 

Evidently, Eqs. (12.671 provide a “Zero Curvature Representation” (ZCR) for the map (I2.17P 
and for the local equations of motion of the discrete quantum evolution system considered in the 
previous subsection. Indeed, as illustrated in Figl3](b), the two-by-two matrices L can be viewed 
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as edge transition matrices for an auxiliary discrete linear problem. Eos. (12.671) then imply that 
products of transition matrices along the paths Q o Q 1 Q 2 and Qo Q\ Q 2 are exactly the same, 
i.e., independent of the choice of path between Q 0 and Q 2 . 

Using the notation (12.401) it is convenient to rewrite the ZCR (12.671) in the form 

Li(A) L+ = L+ Li(A), L 2 (A) = L 2 (A) L~, (2.69) 

where A is arbitrary. 


2.9 Commuting Integrals of Motion 

The phase space of the quantum evolution system (12.601) possesses 4 N degrees of freedom ( 2N 
coordinates and 2N momenta). Let us show that the system has exactly 21V linearly independent 
mutually commuting integrals of motion (IM). Introduce two transfer matrices 

N 

T(A) = Tr (l, (A) L+ L 3 (A) L+ ■ ■ ■ L 2 j V _ 1 (A) L+ A? ) = Tr ]J (l^A) L+) (2.70a) 

n=1 

and 

N 

T(A) =Tr(L r L 2 (A)L 3 - L 4 (A) • • ■ L 2JV (A)) = Tr ]j[ (l 2 " _i L+(A)). (2.70b) 

71 = 1 

It follows from (12.441) and (12.451) that these transfer matrices form a commuting family 

[T(A),T(A')] = [T(A),T(A')] = [T(A),T(A')] =0. (2.71) 

Further, using the zero curvature representation (12.69[) one can easily show that the transfer 
matrices T(A) and T(A) are, in fact, time-independent, i.e., they are unchanged under the 
discrete evolution map (12.581) . 


T(A) = U (T(A)), T(A) = W(T(A)) . 


(2.72) 


Let (Hk, E k , Fk) e Ak be the set of generators belonging to the fc-th algebra in the tensor 
product (12.561) . Define the operators 

V+ = q - 1 E k q -Hi-H 2 -...-H k ^ = q - 1 E k q -H k -H k+ 1 -...-H 2N , 

(2.73) 

V~ = q~ l F k q+Hi+H 2 +...+H k , v fc = q~ l F k q + H k+H k+ 1 +...+H 2N , 

which obey the following relations 


where a, a' 


n n' = q- 2arT ' vf vz, 

vr v; = <? +w v; v;, 

k > £, 

(2.74) 

0 

II 

i> 

[Vfc",v;'] = o, fc + c 


(2.75) 

= +1 and 

[V fc + , V fc - ] = [■v fc + , V-] =(q- q 1 )(q Hk - q~ Hk ) • 


(2.76) 


To within simple factors the transfer matrices (12.701) are IV-th degree polynomials in the 
variable A 2 (or A -2 ), 


N 


N 


T(A) = X N X~ 2n G n , T(A) = A"* £ A 2n G r 


n =0 


n=0 


(2.77) 


15 





















whose coefficients are commuting integrals of motion. 


[G fc , G m ] = [G fc , G m ] = [G fc , G m ] . (2.78) 

From the definition (12.701) it is easy to see that there are only 2N independent coefficients, since 
two pairs of them coincide 

2 N 

G 0 = G 0 = q p + q- F , G n = G n , P = ^F t . (2.79) 

k= 1 

The remaining coefficients can be expressed as “ordered sums” of products of the operators 
(12.731) , which are analogous to the ordered integrals that appear in the context of the continuous 
conformal field theory (33j (the so-called non-local integrals of motion). For instance, the simplest 
non-trivial coefficients read 

N 2n—2 2N 

Gi = - E + <T p+ * a — - « P+1 E v <" v i-i - <T P+1 E v i-i vr} (2-80) 

ra=l t =1 ^=2n 

N 2n—1 23V 

Gr = - 2 {g +p -^ + <r p+i ^ - g p - ] 2 v; V 2 ; - q 1 2 v 2 - v;} (2.81) 

71=1 1 = 1 ^=271 + 1 

Finally note that the transfer matrices (12.701) can be viewed as modified versions of the 
transfer matrix of the six-vertex model |24j . 


2.10 Matrix elements of the quantum R-matrix 


Here we present calculations of the matrix elements of the universal R-matrix for the represen¬ 
tations (12.481) . when the element v is diagonal. Various parts of these calculations previously 
appeared in [34H39] . The resulting expression (12.411) is remarkably simple and essentially coincide 
with that obtained in m- Define new parameters 


q = e 


i7rb 2 


r >=\( b + b ^ 


Im(b 2 ) >0, 


(2.82) 


which are simply related to the value of q < 1. We will also assume |b| < 1. Introduce the 
Heisenberg algebra 


[ x > P] = 


2vr 


Below, it will be convenient to represent the operators entering (12.481) in the form 

v _ g27rb (x+i/3) z _ _ e 2mb(a-/3) 


u = e~ 2nbp 


(2.83) 


(2.84) 


where a and /3 are arbitrary (complex) parameters. Note, that the parameter f3 can be eliminated 
by a trivial redefinition, 


-27Ti/3 E5 


F 


3 + 2771/3 a ^ a+ ^ 


(2.85) 


however, we prefer to retain it for further convenience. Next, consider the coordinate represen¬ 
tations of the Heisenberg algebra (I2.83P in the space of functions /(s), quadratically integrable 


J 


'ZC+OO 


|/(s)| 2 ds < co, 


|Re(c)| ^ |Re(b)| 


( 2 . 86 ) 
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along any shifted real line Re(s) e R in the strip |Im(s)| ^ Re(b), where the operators x and p 
acts as the multiplication and differentiation, respectively, 

x|s> = s|s>, p|s> = |s>. (2.87) 

The corresponding representation of U q (sl(2)), given by (12.481) . (12.841) and (12.871) . is irreducible 
[57] (see also [40]). It will be denoted as Tr a ^. It is spanned by the vectors |s,ct,/ 3), whercH 


v|s,a,/?) = e 27rb ^ s+l/5, |s,a,^), u|s, a, /3) = |s — ib, a, /3), z |s, a, /3) = — e 2mb ( a s,a,/3). 

( 2 . 88 ) 

Further, let 

|si,S 2 > = |si,ai,/3i>® |s 2 ,Q! 2 ,^ 2 > (2.89) 


denotes basic vectors in the tensor product Tf ai ,f} 1 ®^a 2 ,0 2 °f two such representations, where 
for brevity the dependence on the parameters ai j2 and /?i j2 in the LHS of (|2.89l) is not explicitly 
shown. It is useful to introduce additional notations 


_ p27rb(si+i/3i) 

u l c 1 


v '. = e 27rb(s'+i/3i) , 


z — g27rib(a;—/3i) 


i = 1,2. 


(2.90) 


Equations (12.501) and (I2.5ip provide recurrence relations for the matrix elements of the R- 
matrix. For instance, consider the second equation in the left column of (12.511) . namely 


Rv 2 = (—Vi + (v 2 - —vi)ui) R . (2.91) 

V Z2 Z2 ' 

If this equation is “sandwiched” in between the vectors <si,s 2 | and s), s 2 ) it gives 

<si, s 2 | R |si, s 2 > (v 2 - —v^ = (v 2 - —vA (s 1 + ib, s 2 1 R \s[, s 2 ) . (2.92) 

' z 2 / v z 2 / 

In this way four similar equations in (12.511) and (12.531) lead to the following four recurrence 
relations 


<>i + ib,s 2 |R|s / i,4> 

< s i, s 2 | R l s i, s 2> 


Zl 


q 



Z2V' 2 \ 

Zi Vi J 



(s 1 ,s 2 - ib| R|^,4> 

< s i, s 2 | R l s i, s 2 ) 


Z 2 


_ V2_ 

V’l 


1 - 


Z 2 V 2 

qv i 


<Si,s 2 [ R[s( + ib,.s 2 ) 
<Si,s 2 |R|si,s' 2 ) 



<Si,5 2 |R|s / i,s^ - ib> 

< s i, s 2 | R l s i, s 2> 


1 - 


^2^ 2 
Zl Vl 


qz 2 


Vo 


q z i v\ 


(2.93) 

We assume that the matrix elements (s 1? s 2 | Ri 2 |s^, s 2 ) are analytic in the the strip —b ^ 
Imsj,Ims( ^ b. Then the above recurrence relations uniquely define these matrix elements up 
to an inessential normalization factor, 

( s l, ®2 I R-12 I s l 3 s 2 ) = ^/3i-a 2 ( s 2 ~ s l)^ai-/3 2 ( s 2 ~ ai-a 2 ( s 2 — s l)^ 0 1 -/3 2 ( s 2 ~ s l) j (2.94) 

3 Note, that after the redefinition (12.8511 and some other simple equivalence transformations the representation 
7Va,/3 can be reduced to the representation tt 3 defined by Eqs.(2.4) of ref. [33 with s = ia. 
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where 


(2.95) 


, , _ p -i7r/8+i7rs 2 

V a {a) = e l,r ' /8_l7rs tp(\a - s) , P Q (s) = —-:- 

(p[\a — 177 — s 

The function ip(s) is the “non-compact” quantum dilogarithm |36j 

( g ) _ ^ /1 f e~ 2i ^ dw\ 

\4 JiR+io sinh(u;b) sinh(u;/b) w ) 

and the parameter 77 is defined in (12.8211 . Note that essentially the same expression (|2.94D (but 
in a slightly modified form) was also obtained in ref. (37) . 

The non-compact quantum dilogarithm satisfies the quantum pentagon identity [22j 


¥>(pMx) = 7>( x Mp +x)<p(x) , 


(2.97) 


where x, p are the elements of Heisenberg algebra (|2.83l) . 

The functions s V)*(s) and V Q (s) are neither positively defined nor symmetric with respect 
to the substitution s —* —s. However, they are Fourier-dual to each other 


| 

Jr 


-.27T\SX 


V a {s)ds = V a (x) . 


(2.98) 


The factorized form for the kernel of the R-matrix corresponds to a factorized operator 
expression acting in the product of two Heisenberg-Weyl algebra (|2.47|) (but not the full algebra 

U q ( S l( 2 ))), 


R 42 = Vs 1 _ Q2 (x 2 - xi)H ai _ a2 (-pi)l/a 1 _ / 3 2 (p 2 )y ai _ / 3 2 (xi - x 2 )Pi 2 , (2.99) 

where Pi 2 is the permutation operator, 

Pi 2 |si,ai,/3i>® |s 2 ,a 2 ,/3 2 > = |s 2 , aq, /5i> (g) |si, a 2 , fc), (2.100) 

which swaps the coordinates si and s 2 , but keep the parameters ctj and /3{ unchangec0. This 
R-matrix coincides with that for the Dehn twist in quantum Teichmiiller theory [39j. Explicit 
correspondence will be established in Appendix. The elements of the matrix inverse to (12.94|) 
are given by 

^ai- Q2 ( 4 S 2 - S l)^3i-/3 2 ( S 2 “ s ' l) 


<4,41^ 1 |si, s 2 > = 


^1-02 («2 s l)^ai— ( s 2 s l) 


where 


V* a (s) = e^-^VOry + \a - s) . 

This function possesses the property 

lim dsV a+e (s — n)H*_ e (s — v) = 5(u — v ) 


( 2 . 101 ) 

( 2 . 102 ) 

(2.103) 


where a is purely imaginary and a real positive e defines circumventions of the poles. Note, 
that in the case when 77 is real and the spectral parameter a is purely imaginary, the superscript 

4 It should be note that the factorization of the type (12.9411 , (12.991) only happens for a specific choice of the basis 

of representations of U q (sl(2)). It was first observed in [37] for the R-matrix of the chiral Potts model, related to 
cyclic representations of U q (sl(2)). 
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V can be viewed as the complex conjugation. The relation (12.1031) ensures the fact that the 
matrices (12.941) and (12.1011) are indeed mutually inverse. 

The R-matrix (12.941) solves the quantum Yang-Baxter equation!! 


I 


ds[ ds' 2 ds' 3 (s 1: s 2 1 Ri 2 |4> 4) <4; Ail R 13 1*1, 4) <4> 41 r 23 I4> 4) 


r 

Jr 


(2.104) 


ds'i ds 2 ds 3 (s 2 , s 3 | R 23 |4, 4) ( s i> 41 R i 3 14,4) (4 5 41 R i 2 |4> 4)* • 


The statement follows from the Yang-Baxter equation (12.121) for the universal R-matrix (12.141) 
and the irreducibility of the representation of U q (sl( 2 )), given by (|2.48|) and (12.881) . It is worth 
noting that the some solution (I2.94|) can be obtained as a limiting case of the R-matrix of 
the Faddeev-Volkov model mm- Thanks to this connection Ea. (l2.104p can be independently 
verified using the various star-triangle relations, associated with the Faddeev-Volkov model and 
its reductions (see Appendix A for further details). 


2.10.1 Heisenberg evolution operator 

By construction, it is clear that Ea. (|2.60l) describes an Hamiltonian evolution, in the sense 
that there exists an invertible Heisenberg evolution operator U, 

U(0) = UOU' 1 (2.105) 

which reduces the (discrete time) evolution to an internal automorphism of the algebra of ob¬ 
servables. To write this operator explicitly, first remind that we have assumed that the elements 
Zk take only two values, as in (|2.62l) . Next, consider the representation of algebra (I2.56P 

TTq 71"ai,/3i ® ^ct2,^2 ® Tw,/3i ® ^ct2,02 (2.106) 

given by the product of the representations (|2.88l) with alternating values of the parameters a 
and /?. Further, let 

| s> = |si,ai,/3i> ® |s 2 ,Q! 2 ,^ 2 > ® |s3,ai,/3i> ® ■■■ ® |s 2 jv, a 2 , P 2 ) (2.107) 

denotes the corresponding basis vectors. Then the matrix elements of U are given b}l! 

N 

( S | U I S ' ) = <S 2 „_ 1 , S 2n I R I4n+1> 4n) (2.108) 

n =1 

The physical regime corresponds to imaginary values of a*, /J,. In this case R-matrix is unitary 
and therefore the Heisenberg evolution operator is unitary as well. 

To write down the matrix elements of the evolution operator for T ^ 1 steps of the discrete 
time introduce 2 N(T + 1) variables with k = 1,2,... , 2 N and t = 0,1,..., T, such that 

I s j Olt=0) l^4(out) I s , ty\t=T , 1^1, ti &2,ti ••*, &2N,t) (2.109) 

5 The R-matrix (12.941) contains two sets of spectral parameters (ai,/3i) and (q 2 ,/ 32 ), corresponding to the 
two representations appearing in (12.8911 . therefore we write it as (si, s 2 1 R 12 |si, S 2 ). In a similar manner the the 
matrices (sj, s 3 | R 13 |sj, S 3 ) and (s 2 , s 3 | R 23 |s 2 , s 3 ) depend on the third set of spectral parameters belonging to 
the representation Tr a3} p 3 involved in (12.1041) . 

6 The matrix elements of R appearing in this formula are given by (12.9411 with an appropriate substitution of 
the variables Si and s[. The parameters 01 , 2 ,/3i ,2 therein remain unchanged. 
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describes the initial and final states of the system. Then 

r T -1 N T—l 2 N 

Z = (m)'('® I U |s)(out) = I | } } | | { s 2n—l,t i s 2n,t\ | s 2n+l,t+l i s 2n,f)J’ | | ^ | | dSk,t) 

J ( t=0 n=1 t=1 V k=1 

(2.110) 

3 Quasi-classical limit and classical Yang-Baxter map 

3.1 Poisson algebra 

Consider the quasi-classical limit 

q = e i7rb2 — 1, b — 0, (3.1) 

where the parameter b is the same as in (I2.82|) . In this limit the algebra (12.21) becomes a Poisson 
algebra (i.e., a commutative algebra with an associative Poisson bracket { , }, obeying the 
Leibniz rule). Making a substitution 

K - k, E^e, F —* f, (3.2) 

and replacing commutators in (12.21) by Poisson brackets, 

[, ]-27rib 2 {, }, b —> 0, (3.3) 

one obtains the following Poisson algebra, 

V : {k, e} = ke, {k, /} = - kf , {e, /} = k - k~ l . (3.4) 

The central element (12.31) takes the form 

C —> c = ef + k + k- 1 . (3-5) 

The structure of the Hopf algebra in this limit remains intact. In particular, one can define the 
co-multiplication, co-unit and antipode, which, obey the same properties as for the quantum 
algebra U q (sl{ 2)). We will denote them by the symbols A, e and S, respectively. Define the 
co-multiplication A as a map from the Poisson algebra V to its tensor square 

A: V^V®V, (3.6) 

acting as 

A(k) = k(g)k, A(e) = e®/c + l(x)e, A (/) = / ® 1 + k~ l ® / . (3.7) 

Similarly to (12.6|) and (12. 7 1) define the co-unit 6, 

e(k) = l, e(e) = e(/) = 0, (3.8) 

and the antipode S, 

S(k) = k~ 1 , S(e) = -ek-\ S (/) = -k f. (3.9) 

Note that co-multiplication and co-unit are homomorphisms of the Poisson algebra, while the 
antipode is an antihomomorphisms. In particular, this means that 

A ({a, b }) = {A (a), A (6)}, S ({a, b }) = -{S (a), S ( b )} . (3.10) 

The relations (|2.9I) remain exactly the same as in the quantum case, 

5(1) = 1, eo 5 = e, (S®S)oA = A'qS (3.11) 

where A! = a o A is another co-multiplication obtained from (13.71) by interchanging factors in 
the tensor product. 
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3.2 Quasi-classical expansion of the universal R-matrix 

In the limit (13.11) the universal R-matrix (12.141) becomes singular. It is not difficult to show that 


R = (l — e <S> f) 2 exp 


i7rb 2 


(log k (x) log k + Li 2 (e (g) /)) 


(l + 0(b 2 )). (3.12) 


where k, e, f are the generators of the Poisson algebra (13.41) and Li 2 (x) is the Euler dilogarithm 

log(l - t) 


U 2 (x) = - 


f 


■ dt. 


(3.13) 


3.3 Classical Yang-Baxter map 

Even though the quasi-classical limit of the universal R-matrix becomes singular (see (I3.12|) 
above) its adjoint action on the elements of the tensor product of two algebras ( 12 . 21 ) . appearing 
in (12.161) , is well-defined [T7J . Therefore the q —> 1 limit of the quantum Yang-Baxter map (|2.16l) 
is well defined. We will denote it by a special symbol 

1Z = lim 1Z, (3-14) 

to distinguish it from its quantum counterpart. 

Let x = (k , e, /) denotes the set of generators of Poisson algebra V and x\ y2 denotes the 
corresponding sets in the two components of the tensor product V (x) V, 

x\=x®l, x 2 = 1 ® x, x = {k,e,f}. (3.15) 

Similarly to (I2.22p the classical map 1Z is completely defined by is action to two sets of generators 
in the tensor square V ® V, 

7Z: (xi,x 2 ) -» (x^xa), (x[,x' 2 ) = 7Z(x 1 ,x 2 ). (3.16) 

The formulae (12.171) lead to 

r k[ = ki (1 - e\ f 2 k 2 /k 1 ) , 

< e[ = ei k 2 , (3.17a) 

k fi = h/h + h - h K 2 (1 - ei h h/ki)~ 1 , 

and 

r k’ 2 = k 2 (l - ex f 2 k 2 /ki )~ l , 

' e 2 = k\ e 2 + ei - e\ k 2 2 (l - e\ f 2 k 2 /kx) 1 , (3.17b) 

. A = h/ki. 

We would like to stress that this map preserves the tensor product structure of two Poisson 
algebras 

n : V<$)V ^V®V, (3.18) 
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which means that the “primed” quantities {A: ■, e ■, f- }, % = 1,2, have exactly the same Poisson 
brackets (13.41) as the elements {ki,ei,fi}. The inverse map reads 


r kx = k[ (1 - ei/£) \ 

< e\ = e[ {k' 2 )~ l (l - e[ fy) 1 , (3.19a) 

. h = (,/i + f'JK) (1 - ei P 2 ) k' 2 - k\ /' A;' , 

and 

' k 2 = k' 2 (l~ e' x /i), 

< ^2 = (e^ + e i^ 2 ) (l — e i fij/k'i — e'l/ik'i k 2 ) (3.19b) 

^ h = f 2 k[ (l — ei f 2 ) 1 . 

Similarly to the quantum case (see (12.321) ) the direct and inverse maps are related by the antipode 

(ESI) 

o (5 (g) 5) = (5 ® 5) o TZ 1 . (3.20) 


3.4 Properties of the classical Yang-Baxter map 

The quasi-triangular structure of the quantum algebra is also inherited on the classical level. 
Therefore, most of the relations of Sect. El remains literally unchanged in the classical case. 
For completeness we summarize them below. Let 

Xi = {ki,ei, fi}, i = 1,2,3,... (3.21) 

be sets of generating elements of different copies of the Poisson algebra (13.4p . Introduce the 
set-theoretic multiplication 5 which acts on a two sets generators 

5: (xi, x 2 ) x , (3.22) 

and write it as 

x = 5(xi,x 2 ). (3.23) 

Explicitly, one has 

k' = k\ k 2 , e' = e 1 k 2 + e 2 , f = fi + Aq _1 f 2 . (3.24) 

which is essentially the same formulae as (13.71) . As in quantum case, the map (|3.22D is only 
defined on two sets of generating element and the direction of the arrow in (13.22[) is reversed 
with respect to (13.6p . 

Consider the tensor product V ® V (x) V of three Poisson algebras (13.41) and let (x\,x 2 ,x 3 ) 
denote three sets of generators in the corresponding components of the product. Define a 
functional operator 1Z\ 2 

^ 12 (^ 1 , x 2 ,x 3 ) := (R-(xi,x 2 ),x 3 ) (3.25a) 

which acts as (13.161) on the first two sets and does not affect the third one. Similarly, define 
operators 

7^13 = 023 0 7^-12 0 023) Tl 23 = a 22 o TZ 13 o a\ 2 . (3.25b) 
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In the q —» 1 limit Ea. (l2.26p simply reduces to the set-theoretic Yang-Baxter equation for 
the classical maps TZ %] 

TZ 23 ° 77i3 0 ^-12 = 77i2 0 ^13 ° ^-23 • (3.26) 

To formulate the remaining properties of the classical map 77 introduce additional notations, 

5n(x 1 ,x 2 ,x 3 ) := (6(xi,x 2 ),x 3 ) , 5 23 (xi,x 2 ,x 3 ) := (x 1 ,S(x 2 ,x 3 )) . (3.27) 

Then Eqs. (12.11D lead to 

6(x 2 , xi) = 5 o TZ(xi, x 2 ) , (3.28) 

'R(6i 2 (x 1 ,x 2 ,x 3 )) = 6 12 {TZ 23 (7 Zi 3 (x 1 ,x 2 ,x 3 ))) , (3.29) 

'R(5 23 (x 1 ,x 2 ,x 3 )) = 6 23 {TZ 12 ('JZi 3 (x 1 ,x 2 ,x 3 ))) , (3.30) 

where the permutation operator a is defined in (12.101) . Further, from (|2.13l) it follows that 

77 o (e® 1) = (e® 1), 

77 o (1 ® e) = (1 (x) e), (3.31) 

77 o (S (g) S) = (5 ® 5) o TZ 1 , (3.32) 

where the co-unit e and the antipode S are defined in (13.91) . We would like to stress that the 
above relations (13.261) . (13.2811 (13.32 jl are not specific to the Poisson algebra (13.41) and must hold 
for any (quadratic) Poisson algebra which arise in the quasi-classical limit of a quasi-triangular 
Hopf algebra. In the case under consideration these relations could be verified using the explicit 
form of the map 77 and its inverse given in (12.171) and (12.181) . Note that, contrary to the 
quantum case, the above relations only involve ordinary substitutions for a set of commuting 
variables, since the Poisson algebra V is a commutative algebra. The fact, that the map (13.181) 
is a “canonical transformation” (in the sense of Hamiltonian dynamics), preserving the Poisson 
brackets of the generating elements, does not play any role for the validity of (I3.26l) - (13.32l) . 


3.5 r-matrix form of the Poisson algebra 


As noted before, the universal R-matrix (|3.12j) as well as the Yang-Baxter equations (12.121) and 
(12.341) becomes singular in the quasiclassical limit (13.11) . (13.21) . We will return to these singular 
equations in more details in Sect l3.1(Tl 

Here we consider the remaining equations in Sect 12.51 which all have a smooth limit q —> 1. 
The classical analog of the L-operators (|2.37|) and ()2.40|) reads 


//c2 kz f\ 

(o * 



and 

£{X) = X£ + - A -1 £ - . 

For the R-matrix ()2.39|) one obtains 



(3.33) 


(3.34) 


R ± = 1 ± 


i7rb _ 


2 i7rb 2 r~ 


0( b 4 ), 


0 , 


(3.35) 
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where 


/ ° ; \ 
i 


i _i 

I 2 

V : 0/ 

Similarly to (12.41 p and (12.431) define 


/ o 


r = 


i 

i 

1 i 

2 i 
_L 


-1 


I 1 


0 / 


(3.36) 


r(A) = Ar + — A 1 r , r(A) = r(A) P (3.37) 

where P is the permutation operator and r(A) is the classical r-matrix satisfying the classical 
Yang-Baxter equation |42| 

[n 2 (A), n 3 (A/r)] + [ri 2 (A), r 23 (g)] + [ri 3 (A/u), r 23 (/r)] = 0. (3.38) 

The equation (12.451) then reduces to 


{£(A) (g) £(g)} = [r(X/ii), £(\)®£(^)] , (3.39) 


where the symbol { (x) } denotes Poisson bracket of two factor of the tensor product. 


3.6 Heisenberg-Weyl realisation (canonical variables) 

The three-dimensional Poisson algebra (13.41) has one central element (13.51) . Therefore it must 
essentially reduce to the classical Heisenberg-Weyl algebra, defined by the bracket 


W : {u,v} = uv. (3.40) 

The required coordinate transformation (cf. ()2.48|) 1 is 

k = u , e = v (z — u) , f = u -1 (l — z~ l u~ l ) . (3.41) 

The variable z parametrises the central element (13.51) . 

c = z + z _1 , (3.42) 

and has vanishing Poisson brackets with u and v. In the new variables the map (13.171) reduces 
to the classical analog of (1 2.51 [) 

u[=uig c i, u' 2 = u 2 /g c i, 

< = V 1 V 2 U 2 (ui u 2 + (u 2 — V 1 /Z 2 )) 1 V 2 = z\Vi/z 2 + {v 2 - Vi/z 2 )ui, (3.43) 

z[ = z 1 , z 2 = z 2 


where 


9d = 1 - 


Vl <yZ\ 


ui){u 2 - 1 / 2 : 2 ) 

U\V 2 


(3.44) 
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Note that the parameters z\ 2 remain unchanged under the map. Similarly for the inverse map 
one obtains 


“i = u i/9d, 


u 2 = u 2 g'di 


Vl 


v[v ' 2 (ziv[/z 2 + (v ' 2 - Ziv[) U 2 ) 


-1 


V 2 = v[ + (v 2 ~ Z\V\)/vl\ . 


(3.45) 


where 

j n v[(zi-u' 1 )(u' 2 -1/z2) 
9d ~ 1 -J—, - 


(3.46) 


3.7 Discrete classical evolution system 


The quantum evolution system considered in Section [2 .71 admits a straightforward classical limit. 
Remind, that the classical Yang-Baxter map (13.161) acts on the tensor square of the Poisson 
algebra (13.41) . 

77: (Vi® V 2 ) ~ (V\ (x) V' 2 ) = n(Vi ® V 2 ) • (3.47) 

Correspondingly, the algebra of observables (12.561) becomes a tensor power Poisson algebra (13.41) . 

O c i = V\ ®V 2 ® ■■■ ®V2N-i ®V 2 n, TV ^ 1. (3.48) 


Similarly to the quantum case the factors of this product are assigned to edges of one horizontal 
saw, corresponding to a fixed moment of the discrete time, as shown in Fig. [4] (one just need to 
replace all .47 s there by V’s to get the classical picture). The classical evolution map is defined 
as ^ ^ ^ 

Mel = S o (lZi 2 o 7734 o • • • o 77 ( 2 Tv-i), 2 iv) > (3.49) 

where 77 ij = ay,- o77jj is the map (13.471) , acting on the 7-th and j-th factors of the product (13.481) , 
followed by the permutation operator oy,-. The operator S cyclically shifts (to the right) the 
factors in the product (13.48p . similarly to (I2.59p . The one-step time evolution is defined as 

O r c i = U c i(O c i) = V[ (x) V 2 ® • ■ ■ ® V 2N _i ® V 2N . (3.50) 

For the realisation (I3.4ip each edge will be assigned a triple of dynamical variables {ui, Vi,Zi}, 
i = 1, 2,... , 27V. As in the quantum case we assume the homogeneous arrangement (12.621) . Then 
all z’s stay unchanged under the time evolution. The remaining variables possess the canonical 
Poisson brackets 

{\ogUi ,\ogVj } = <%, i, j = 1,2,... ,27V. (3.51) 

which are preserved by the time evolution. 

Taking classical limit in (12.631) and introducing discrete time variable t , one immediately 
obtains the Hamiltonian equations of motion 

^2n+l,t+l = 'U j 2n—l,t9n,ti 


where 


1 


^2n+l,t+l = 

,V2 n,t 

U 2n,t+l = 'U'2n,t/9n,t j 
Zl 


1 


1 


1 


Z2 


^2n—l,t Z2V 2 n,t J ^2 n,t 


1 


V2n,t+1 V2n—l,t T I V2n,t '^2n—l ) t j ^2n—l,t ; 


Z 2 


9n,t — 


l _ V2n-l,t( z l ~ U2n-l,t)(u2n,t ~ V^) 
^2n—l,t^2 n,t 


(3.52) 


(3.53) 
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3.8 Zero curvature representation 

The quantum zero curvature representation (|2.67|) admits a straightforward quasi-classical limit. 
Let if an d £2 denote the classical ^-operators (13.331) with elements belonging respectively to 
the first and second Poisson algebra of the tensor product V\ (x) V 2 - Next, let if and if are 
defined in a similar way, but with elements from the “transformed” algebras V[ and arising 
from the map (13.47|) . Then Eas. (l2.67p lead to the relations 


it it = io T 


+ t;- 


1 ) 


If it =io i 


1 5 


if if = in i 


l > 


(3.54) 


which can be easily verified using explicit formulae for the classical map (|3.17l) . Conversely, 
the equations (13.541) can be used as an alternative definition of the classical Yang-Baxter map 
independent of the notion of the universal R-matrix for the quantum algebra. The geometric 
interpretation of the relations (|3.54|) as the zero curvature representation remains essentially the 
same as in the quantum case, see Figj3](b) and the relevant discussion at the end of Sect l2.8l 
Finally, note that Eas. (|3.54l) can be conveniently rewritten as 

4(A) it =it 4(A), it 4(A) = 4(A) if , (3.55) 


where A is arbitrary, similarly to (12.69[) . 


3.9 Involutive Integrals of Motion 

The following considerations essentially repeat those for the quantum case in Sect J2.9l The phase 
space of the classical evolution system (13.501) possesses 4 N degrees of freedom (2N coordinates 
and 2 N momenta). Let us show that the system has exactly 2 N integrals of motion (IM), which 
are in involution to each other with respect to the Poisson bracket (13.41) . The classical analog of 
the transfer matrix in the quantum case is the trace of the monodromy matrix for an auxiliary 
linear problem. Similarly to (|2.70D introduce two such traces 

N 

i(A) = TV (i!(X)i+ 4(A) it • • • 4jv-i(A)^+v) = TV (i 2n -i{\)i£ n ) , (3.56a) 

71 = 1 


t(A) = Tr (if 4(A) if 4(A) ■ ■ • if N _i 4v(A)) = TV fj fc n -i4n(A)) • (3.56b) 

n =1 

It follows from (13.391) that these quantities form an involutive family 

{t(A),t(A')} = {i(A),t(A / )} = {i(A),t(A')} = 0, (3.57) 

for arbitrary values of A and Ah 

Further, using the zero curvature representation (12.691) one can easily show that t(A) and 
t(A) are, indeed, time-independent, i.e., they are unchanged under the discrete evolution map 

M, 

t{ A) = Ud (t(A)) , t( A) = U d (t(A)) . (3.58) 

Similarly to (12.771) in the quantum case, the classical Integrals of Motion can be defined as 
coefficients of expansions of the polynomials i(A) and t(A) in the variable A 2 . 
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3.10 The Lagrangian equation of motion and the action 

In view of the quasi-classical correspondence, considered in the section, the Lagrangian function 
and the action of the classical evolution system (13.521) is determined quasi-classical asymptotics of 
the universal R-matrix (13.121) . More specifically, consider the expression for its matrix elements 
given by Eq. (12.941) . It involves the spectral parameters on, 2 , /?i ,2 and the spin variables sq 2 and 
s) 2 - For a consistent quasi-classical limit these variables should become infinite, provided that 
the quantities 


a,i = — 27ribaj, bi = — 27rib/3j, <7j = 27rbsj, a[ = 27rbs(, i = 1,2, b —> 0, (3.59) 

remain finite for b — * 0. Introduce new functions 

2 2 

Aa(cx) = -y - Li 2 (-e- CT -“), A a (cr) = y + Li 2 (e~ CT -“) , (3.60) 

where the Euler dilogarithm Li 2 (x) defined in ()3.131) . It is useful to note that 

d -^l = - log(e CT + e-“), = log(e CT - e” 0 ) . (3.61) 

With these notations one can show that 

(s 1 ,s 2 \R\s' 1 ,s' 2 ) = exp + 0(1)| as b -> 0 , (3.62) 

where 


C — Ab 1 _ a2 (fJ 2 — a 1 ) + A ai _b 2 ((J 2 — (j'i) + A ai _ a2 (cr 2 — aq) + —6 2 (cr 2 — <?'{) ■ (3.63) 

This formula easily follows from the asymptotics of the functions (12.951) for b —* 0, 


27rb 


(7 

27rb 


= exp 


27rb 2 


Aa(<7") + 0(1) 


-rv ( <J 


V 


27rb 


= exp 


27T b 2 


Aa( 0 ') + 0(1) 


(3.64) 

Eqs (|3.52l) are the (first order) Hamiltonian equations of motion for our dynamical system. 
To obtain the corresponding (second order) Lagrangian equations of motion one may consider the 
quasi-classical asymptotics of the matrix elements of the evolution operator (12.1101) . Substituting 
(13.621) therein and calculating the integral by the saddle point method one gets 

1 °gZ = ^A{{a}W)+°( 1 ) (3.65) 


where M({(j}^ c ^) is the classical action of the system 

N T—l 

AM) -EE £(02ra-l,t) a 2n,ti cr 2n+l,t+l, &2n,t+ 1) (3.66) 

n=l t =0 


calculated on the stationary spin configuration {a}^ defined as 

dA({a}) 


8(TkA 


[cr} = {cr}C c ^) 


= 0. 


k = 1,..., 2N, t = 1,..., T — 1. 


(3.67) 


The boundary spins {&k,o}k=i and Wk.T}'k=\ are kept fixed 


2N 
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The same result can be obtained by the standard method of Hamiltonian mechanics. Remind 
that the classical map (13.431) involves the set of six independent variables 1 * 1 , 2 , ^ 1 , 2 , ^ 1 , 2 - They 
are simply related to the variables in (13.591) . 


_ cfii a i 


Zi = —e 


Vi = e ai bi , 



(3.68) 


Next, the map (13.431) is a canonical transformation, preserving the Poisson bracket 


{log u i} log?;,,} = <%, i,j = 1,2 


Its generation function L is defined by the equations 

dC= ^ (logu'idlogv'i-logUidlogVi) . 

i= 1,2 


or 


SC 

d log Vi 


Rewrite (13.43 j) in the form 


= log V>i , 


dC 

d log v[ 


log u'i, 


i = 1,2. 


u 1 = z 1 
u[ = z 1 


1 - z 2 v' 2 /z 1 v 1 

1 — Z 2 V 2 /VI 
1 - v'Jzlv^ 

1 - v 2 /v[ 


u 2 = 

u 2 = 



1 - Z2V2/V1 

1 - v 2 /v[ 

1 - z 2 v' 2 /ziVi 
1 - v' 2 /z 1 v[ 


(3.69) 


(3.70) 


(3.71) 


(3.72) 


Integrating (|3.70l) and using the variables (13.681) . one obtains the same Lagrangian (13.631) . 



Figure 5: Arrangement of classical variables v 2n -i,t and v 2n ,t- 


3.11 Discrete Liouville Equations 

It is instructive to write down the Lagrangian equations of motion (I3.67P explicitly. The 
corresponding Hamiltonian equations (13.521) are written in terms of the variables Uk,t-> v k,ti 
k = 1,...,2 N, where the variables v^t (playing the role of coordinates) are related to ak,t 
in (I3.66P as follows 

V2 n-M = e- 6l+<r2 "- 1 ’*, V 2 n,t = e- b2+r72 ^ , n = 1, ... ,N, t = 0,...,T. (3.73) 
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They arrangement around an elementary quadrilateral is shown in FigEJ Rewriting (13.431) in 
the form 13.721 one obtains 


U2n-l,t 


zi- 


1 _ Zi V2n,t+1 
%2 v 2n-\,t 


1-^2 


^2 n,t 
v 2n-l,t 


U2n,t = 


1 ^2 n,t 

, 1-22 - 

1 ^2n-l,t 


22 1 _ 


V2 n,t 


t>2n+l,t+l 


U2n+l,t+l = 2l 


1 ^2n,t+l 

2l V2n+l,t+l 
V2n,t 


^2n+l,t+l 


U2n,t+1 


j _ 22 V2n,t+1 
1 Zi V2n-l,t 

Z2 -y 1 '^2n,t+l 

2l V2n+l,t+l 


(3.74) 


To obtain the Lagrangian equations of motion (13.671) one needs to exclude the variables u^i 
(they are regarded as momentum variables). This is achieved by equating u’s defined by the 
first line of (13.741) to those from the second line of (|3.74l) (with an appropriate shift of indices). 
Consider a geometric structure of the resulting equations. Let v = v^t be one of the v -variables 




Figure 6: Arrangement of variables for the Lagrangian equations of motion, 
for some fixed values of k,l. Define four adjacent variables 

Di = Vk+ijt+ii Vd = 'ttfc—i,t— 15 ty = Vk—i,ti ty = Ufe+gt, v = , (3.75) 

shifted one site “up”, “down”, “left” or “right” with respect to v, as shown in FigEJ The 
equations (13.671) has different form depending on whether the spatial coordinate of the central 
site (k, t ) is odd or even. For an odd k = 2n — 1 (left part of Figj6|) one obtains 


and for an even k 


(i _ (\ _ Z2Vr 

V z±vJ V v 
2 n (right part of FigEJ), 




22 v \ 
21 v d ) 



(3.76a) 


(3.76b) 
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3.11.1 General solution of the equation of motion 

Below we present a general solution of both the Lagrangian and Hamiltonian equations of motion. 
It will be convenient to introduce the “light cone” coordinates: 


x = (n,t — n), ei = (1,0), e 2 = (0,1). 


Conversely, 

n = X 1 , t = X 1 + X2, X=(xi,X 2 ). 

Moreover, we will use the variables ui j2 (x) and ug 2 (x) instead of Uk,t and Vk.t, 

Ui(x) = U2n-l,t, U 2 {x.) = U2 n ,ti «l(x) = V 2n -1 ,t, « 2 (x) = V 2ri)t • 


(3.77) 

(3.78) 


(3.79) 


The Hamiltonian equations of motion in these notations are easily obtained from (|3.72l) by 
substituting the variables {vi,v 2 ,v[,V 2 } therein with (ui(x),u 2 (x),tq(x + ei),u 2 (x + ei)} and 
similarly for {ui, u 2 , u'i, u' 2 }- Their arrangement on lattice is shown in Fig. [3 (cf. Fig. [5]). The 


T x+ei 


r x+ei+e2 



Figure 7: Light cone coordinates, x = x\e\ + x 2 & 2 - 
resulting equations imply a simple relation 

ui(x)u 2 (x) = iti(x + ei)tt 2 (x + e 2 ) . (3.80) 

which ensures the existence of the r-function, dehned by a system of two first order difference 
equation 

ttl(x) = - , u 2 (x) = - . 

'T’x ^x+ei 

With this substitution the system (13.741) can be rewritten as 

ui(x + ei) = zit x - r x+e2 

«l(x) ■SlLx+ei 7x+ei+e2 

^2(x + e 2 ) _ ^'2' r x+e2 T x+ei+e2 

v 2 (x) z 2 t x - r x+ei 

Wl(x) _ ^ 7 ~x+ei’bc+e 2 7 ~x7x+ei+e2 

U 2 (x) (ziT x - T x+e2 )(z 2 T x - r x+ei ) ’ 


(3.81) 


(3.82) 
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Here we have only three equations, since (13.801) is satisfied automatically by virtue of (13.811) 
Integrating the first two equations, one obtains 

Px 2 


ui(x) = 


, v 2 (x) = 


Z2T x T x | ei 
%2 


ZlT x T x _|_0 2 

where x = (£ 1 , 0 : 2 ) are the light cone coordinates, 

x = aqei + £ 2 e 2 , x\ = n , x 2 = t — n , 


(3.83) 

(3.84) 


and a Xl , /3 X2 are arbitrary constants of integration. The third equation in (I3.82|) reduces to a 
single equation for the r-function, 


’^"x+ei^”x+e2 7x7x+ei+e2 ^x\ 


Px 2 ■ 


(3.85) 


This is the inhomogeneous discrete Liouville equation, involving two arbitrary functions a Xl and 
j3 X2 . The corresponding geometric arrangement of the T-functions is shown in FiglT] This is a 
generalization of the homogeneous discrete Liouville equation (with a Xl = 1, j3 X2 = 1), previously 
studied in [T9II22] . It appears that the inhomogeneity does not bring notable complications. 
Indeed, a general solution to Eg. (13.851) is given by the formula 


1 + fx\9x 2 
4 > x\'Yx 2 


(3.86) 


where <j> Xl , 7 X2 are arbitrary functions, while f Xl and g X2 are determined by the following first 
order difference equations 


fx 1 + 1 fx 1 — G^l^Xl^Xl + l ) 5x2 + 1 9x2 — /^X2'7x27x2+1 ■ (3.87) 

To summarize the formulae ((3.81 j> , (13.831) , (13.861) and (|3.87[> , together with the definitions (13.731) 
and (13.791) . provide the most general solution to both the Lagrangian equations (|3.67l) . (13.761) 
as well as the Hamiltonian equations of motion (|3.74j) . The solution contains 41V arbitrary 
constants a Xl , {3 X2 , (j) Xl , j X2 (x\,x 2 = whose number precisely coincides with the 

dimension of the phase space of our classical discrete evolution system (I3.50p . 


4 Conclusion 

In this paper we have shown that the theory of integrable maps (and more specifically the so- 
called Yang-Baxter maps) can be naturally included as a part of the modern theory of integrable 
system, based on the theory of quantum groups, the Quantum Inverse Problem Method and 
their classical counterparts. The considerations apply to both quantum and classical theories 
(the latter arise as the quasiclassical limit of the quantum case). One of the main advantages of 
our algebraic approach is that the resulting discrete integrable evolution systems automatically 
possess meaningful Hamiltonian structures. In addition, the problem of quantization of the 
Yang-Baxter maps appears to be solved from the very beginning. 

We have illustrated the entire scheme on the example on the quantized Lie algebra U q (sl( 2)) 
(and its Poisson algebra limit for the classical case). Completely parallel presentations of the 
quantum and classical cases are given in Sect. 2 and Sect. 3 respectively. In particular, a general 
solution the classical discrete Liouville equations arising in this context (both in the Lagrangian 
and Hamiltonian forms) is given in Sect. 3.11. 

Although the above considerations were mainly restricted to the case of U q (sl{ 2)), our ap¬ 
proach is rather general and can be applied to any quantized Lie (super) algebra as well as their 
affine extensions. We hope to consider some of these problems in future publications. 


31 































Acknowledgements 

We thank Ludwig Faddeev, Rinat Kashaev and Joerg Teschner for important comments and 
Nicolai Reshetikhin for alerting us about Sklyanin’s paper m We also thank Sergei Khoroshkin 
and Zengo Tsuboi for numerous fruitful discussions and collaboration on related projects. 


Appendix A. Additional properties of the R-matrix 


Here we present some additional relations connected with the R-matrix (I2.94|l . The non-compact 
quantum dilogarithm (|2.96l) (sometimes called the Barnes double-sine function) is a unique 
solution of the functional equation 


log ^(2 - '\) - log v(z + i^) = log(l + e 2nbz ) 
provided that ip(z) is analytic in the strip 

— Re(b) ^ Im(z) ^ Re(b) . 

It satisfies the relations 

M-z) - 1 , V ( z y , -L. , 

Define also another function 

$ (z)-expf 1 f e ~ 2 ' ZW dw \ 

\8 J K+i0 sinh(u;b) smh(u)/b) cosh(2rjw) w ) 


(A.l) 


(A-2) 


(A.3) 


(A.4) 


which possesses the properties 

$(*)$(-*) = e^ 2 / 2 -^^- 8 ^ 2 )/ 12 , ^(^f)| JB _ J> _ 00 —► 1 . (A.5) 

The Boltzmann weights for the Faddeev-Volkov model (which is related to the quantized 
affine algebra U q (sl2)) are defined as [38] . 

W a (s) = , W a (s) = W v - Q (s) , (A.6) 

F a ip(s - 1 a) 

where 

F a = e i7r “ 2 + i7r ( 1 - 8? ? 2 )/ 24 $(2ia) . (A.7) 

It is useful to note that W a (s) = W a {—s). Consider also a simple limit, 

s —> s — K , ia —► i a — K , K —> 00 (A.8) 

In this case, from (1A.3I) and (IA.5I) it follows that 

W a (s) — e~' mK2+27T ' ,Ks V a (s ) , W a {s) — e mK2 - 27TiKs V a (s) , (A.9) 

where V a (s) and V a (s) are defined in (I2.95|) . 

The weights W a (s ) and W a (s) satisfy the following star-triangle relation [38li4lj . 

f daW a (a-a)W a+l3 {c-a)Wp{b-a) = Wp(a - c)W a+ p(a - b)W a (c - b) (A.10) 

Jr 
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which admits several interesting limits. For example, if 

a —» a — K , c —* c — K , \a —> \a — K, K —> oo , (A.11) 

then the singular terms in (IA.10I) cancel out and the relation reduces to 

f daV a (a-a)V a+ p(c-a)Wp(b-a) = Wp(a-c)V a+p (a-b)V a (c-b) (A.12) 

Jr 

A similar limit, 

a—> a + K , c —> c + AT , ia —► i a — K, AT —» oo , (A.13) 

yields 


f 

Jr 


da V a (a - a)V a+/3 (a - c)Wp(b - a) = W /3 (a - c)V a +p(b - a)V a (b - c) 


(A.14) 


These relations can be recursively used to prove the Yang-Baxter equation (12. 104jl given in the 
main text. 

It is interesting to study the classical limit of the above relations (1A.12I) and (1A.14I) . The 
classical limits of functions V and V is given by (13.6411 . Similar limits of W and W are given by 


Wja_ 

27Tb 



exp 


27rb 2 


A a (&) + 0(1) 


Wja_ 

2nb 



exp 


i — 


27rb 2 


A a (a)+0(1) 


(A.15) 


where 


Then, for instance, 
by the relation 


dA e a + e a dA e a ~ e a 

~da ~ ° g 1 + e a+a ’ ~da ~ ° S 1 - e CT -“ ' 

the stationary point of the integral (1A. 141) in the limit b 

3 

— (\ a (a - a) + \ a+ p(a - c) + Ap(a -b)) = 0 , 


(A.16) 

0 is determined 
(A.17) 


which is equivalent to 

log( e -“ - e~ a ) - log(e CT-c + e~ a ^) + log = 0 • (A.18) 

The last relation has the structure of the Adler-Bobenko-Suris m tree-legs equation 

^(xo, yo; ao) - V’fao, 2/i5 «i) = <p(x 0 ,xr,a 0 , aq) (A.19) 

corresponding to the if 3 system in their classification. 


Appendix B. Relation to R-matrix for quantum Teichmiiller the¬ 
ory 

The R-matrix (12.991) coincides with the inverse of R-matrix (44) from [[39]. Details are the 
following. Moving quadratic exponents in (|2.99l) to the right, one could obtain the left hand 
side of the following identity: 

j e -i7r(xi-x 2 ) 2 e -i7r(p2+p2) e -i 7r (x 1 -x 2 ) 2 p i;i = e -27ripip 2 _ (B.l) 


33 

























This identity can be easily proven by consideration of a kernel or of a normal symbol of the left 
and right hand sides and performing the Gauss integration. 

Using the Pentagon identity (|2.97p twice in (12.991) . one can then bring R 12 to the form 

R12 = F 12 y?(i(/3i - a 2 ) + xi - x 2 )F^e- MpiP2 , (B.2) 

where 

F i 2 = </?('(«! - P 1 ) + PiM '(«2 - P 2 ) ~ P 2 ) ■ (B.3) 

This form of R 12 literally coincides with R -1 , formula (44) from [39] . after indentifying corre¬ 
sponding variables. The authors are grateful to R. M. Kashaev for this observation. 
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